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'"The ALEPH data on hadronic r-decay is throughly analysed in the framework of QCD. The per- 
tiff^ative calculations are performed in 1-4-loop approximation. The analytical properties of the polar- 
ization operators are used in the whole complex q 2 plane. It is shown that the QCD prediction for R r 
ag^es with the measured value R T not only for conventional Af nv = (618 ± 29) MeV but as well as for 
A^" = (1666 ± 7) MeV. The polarization operator calculated using the renormgroup has nonphysical 
ctQT\[— A3, 0]. If A3 = A^ onv , the contribution of only physical cut is deficient in the explanation of the 
A^CPH experiment. If A3 = A^ ew the contribution of nonphysical cut is very small and only the physical 
ctlCexplains the ALEPH experiment. The new sum rules which follow only from analytical properties 
o^3larization operators are obtained. Basing on the sum rules obtained, it is shown that there is an es- 
sential disagreement between QCD perturbation theory and the r-lepton hadronic decay experiment at 
conventional value A3. In the evolution upwards to larger energies the matching of r(q 2 ) (Eq.(12)) at the 



m^ses J/ip, T and 2m t was performed. The obtained value a s (— m 2 ) = 0.141 ± 0.004 (at A3 = A^ ew ) 
difi&rs essentially from conventional value, but the calculation of the values R(s) = —, j -^ hadrons ) 



<r(e+e — »/^ + ^ ) ' 

■R;k = ; F r(z^ieptons) ' a s(— 3 GeV 2 ), a s (— 2.5 GeV 2 ) does not contradict the experiments. 

■ X ■ 

l£ INTRODUCTION 

The purpose of this work is to combine the analyticity requirements of QCD polarization operators 
with the renormalization group. This work is the continuation of works [1-3]. In the work [1] analytical 
properties of polarization operators were used to improve perturbation theory in QCD. In the works 
[2,3] the high precision data on hadronic r-decay obtained by the ALEPH [4], OPAL [5] and CLEO 
[6] Collaborations were analyzed in the framework of QCD. The analyticity requirements of the QCD 
polarization operators follow from the microcausality and the unitarity, therefore we have no doubts 
about them. On the other hand, the calculation according to renormgroup leads to appearance of 
nonphysical singularities. So, the one-loop calculation gives a nonphysical pole, while in the calculation 
in a larger number of loops the pole disappears, but a nonphysical cut appears, [— A§, 0]. As will be 
shown, there are only two values of A 3 , such that theoretical predictions of QCD for R T y+A (formulae 
(23), (24) agree with the experimens [4-6]. These values are the following: one conventional value 
A conv = ( 61g ± 29) MeV and the other value of A 3 is A% ew = (1666 ± 7) MeV. Only in these values 
of A3 the predictions of QCD are consistent with the experiments [4-6]. As far as I know, the value 
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A 3 = A% ew was not considered before now 1 If one simply puts out the nonphysical cut and leaves the 
conventional value A^ onv , then the discrepancy between the theory and experiment will arise. As will 
be shown, if instead of the conventional value Af nv one chooses the value A^ ew = (1565 ± 193) MeV 
then only the physical cut contribution is enough to explain the experiment of the hadronic r-decay. 2 . 
It is convenient to introduce the Adler function (11-13) instead of the polarization operator. The Adler 
function is an analytical function of q 2 in the whole complex q 2 plane with a cut along the positive q 2 
semiaxes. We will use the renormgroup only for negative q 2 , where the value a s (q 2 ) is real and positive. 
For other q 2 the value a s (q 2 ) becomes complex and is obtained by analytical continuation. 
The plan of the paper is the following. 

In Section 2 the formulae obtained in paper [3] are transformed to suitable for this paper form. 

In Section 3 the values A 3 are found such that R T y + A = 3.475 ± 0.022 (formula (24)). 

In Section 4 we obtain new sum rules for polarization operator which follow only from analytical 
properties of the polarization operator. These sum rules imply that there is an essential discrepancy 
between perturbation theory in QCD and the experiment in hadronic r decay at conventional value of 
A 3 . The power corrections and instantons cannot eliminate this discrepancy. 

Section 5 suggests the method of resolving these discrepancies. At A 3 = (1565 ± 193) MeV the 
nonphysical cut gives no contribution into R T y + A and the physical cut gives the experimentally observed 
value R T y + A- The previously derived sum rules make no sense if A 3 is as large. 

In Section 6 we go over to larger energies. In the matching procedure we require continuity of r(s) 
(Eq.(12)) [1] at masses J/ip, T and 2mt when going over from rif to rif + 1 flavours. The number of 
flavours on the cut is a good quantum number. At every point off the cut all flavours give a contribution. 
This follows from the dispersion relation for the Adler function. The continuity requirement off the 
cut when changing the number of flavours violates analytical properties of the polarization operator. 
Section 6 presents the results of the calculations in 1-4 loops for estimation of the precision of the 
calculations. In Sections 7-9 we compare the theory with experiment. In Section 7 the prediction of 
the function R(s) is compared with experiments. The calculated values of the function R(s) are in a 
very good agreement with the experiment (Tables 5,6) at 2 < ^/s < 47.6GeV except for the resonance 
region. 

In Section 8 we compare the calculated values of a s (—3GeV 2 ) and a s (—2.5GeV 2 ) with the values 
a s (—3GeV 2 ) and a s (— 2.5GeV 2 ) obtained from the Gross-Llewellyn-Smith sum rule [7] and the Bjorken 
sum rule [8]. The results of the calculations are in agreement with the values obtained from the 
experiment using these sum rules. 

Section 9 presents the calculation of Ri = ^(z^itptonl) ■ ^ ne obtained value Ri does not contradict 
the experiment. 

Section 10 is devoted to discussion on the analiticity a s (q 2 ). It is shown that the statement that 
a s (0) = 4:ii/Po = 1.396 is valid only for one-loop calculations. 



^^Onc can see this also from Fig. 2 of ref.[3], the line obtained in conventional approach crosses the experimentally 
allowed strip of R T y + A at two values of A 3 equal to A^ onv and A% ew (in rcf.[3] the parameter a s (— m^) related to A3 was 
used). In ref.3 the new value A^ ew was not considered. 

2 The errors here and the following formulae are due only to the error in the measurement of the value R r y + A ( e q-(24)) 
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2 INITIAL FORMULAE 

In this section the formulae obtained in the paper [3] are transformed to suitable for this paper form. 
We will consider three loop approximation thoroughly. Polarization operators of hadronic currents are 
defined by the formula 

KM) = * J e tqx (0 | T J^x) J+(0) | 0}d*x = 

= (<wz* - g»vq 2 M\ q 2 ) + wJlWtf) (1) 

where 

J = V, A; V^ = u-y^d, = u^^d 

Imaginary parts II j\ II j ^ are connected with the measured, so called spectral functions 
vi(s), di(s), do(s) by the formulae 

«i(s)/ai(s) = 2tt ImU { yj A (s + iO) 

a (s) = 2tt ImUf (s + iO) (2) 

Functions II^L are analytical functions of g 2 with the cuts [4m 2 , oo] for 11^ and [9m 2 , oo] for 

IT?, a (s) = 2n 2 f 2 J{s - m 2 ), f n = 130.7MeV. 

To get QCD predictions let us use the renormalization group equation in 3-loop approximation [9-10] 

q 2 ^ = -(3t f) a 2 {l + bf j) a + bt }) a 2 ) (3) 

where 

a ( a = ^, flS"'» = ll-|»„ = 51 - f „„ 

ri g / 27 / 1 ^ n /) 2/5o™ 

Here n/ is the number of flavours. 

Let us consider for the moment rif — 3 and omit the mark rif. Find singularities of a(q 2 ). Integrate 
equation (3) [3] 

(3 ln ^ = - / 2 , a 2 . (5) 

«0 2 ) 

Denote the value q 2 at which a(q 2 ) = oo as — A§. 3 Then we get instead of eq.5: 

(2 \ 00 7 
= / a = /( Q ) (6) 

A§ / J a 2 (l + 6ia + b 2 a 2 ) w v 7 

a(g 2 ) 

According to the known value a(— m 2 ) A3 is determined by the formula 
3 This is the definition of A3. 
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Ha(-mj)) 

A3 = m 2 e Co 



The integral in the formula (6) is taken and the answer is written as 



f(a) = — h 61 Ina . 

a ^2 _ Ah2 [x\ 



-2 £n(a — x\) — — ln(a — x 2 ) 



(7) 



(8) 



where 



Xl,2 



-61 ± y^ 2 - 46 2 
2&o" 



(8a) 



At a s (-m?) = 0.355 [3] we obtain A 2 = 0.39AGeV 2 . 

The expansion of the function f(a) in the Taylor series at large a over 1 /a is of the form 

f(fl) = il_il + !!iMl + o ( l) 

/W 3b 2 a3 46|o 4+ 561 a5 + U V 
It follows from eqs.(6-9) that the singularity a(q 2 ) at g 2 — > —A3 has the form [3] 



(9) 



a(g 2 



2A 2 



1/3 



(10) 



V 3/5 2 (g 2 + A|); 

Since for massless quarks the contributions of V and A coincide, we will omit in all formulae the mark 
J. Introduce the Adler function 



dq 2 



It is convenient to write for three flavours 



D(q 2 ) = 3(1 + d(q 2 )), R(q 2 ) = 3(1 + r(q 2 )), U(q 2 ) = 3(-Zn(-g 2 /V) + p(q 2 )) 

r(q 2 ) = -Imp(q 2 ) = ^-[p(q 2 + t0) - p(q 2 - i0)} 
it 2m 



(12) 



In 3-loop approximation for MS renormalization scheme function d(q 2 ) for negative q 2 is written as 
[11]: 

d{q 2 ) = 4a(g 2 )(l + Ad^a^q 2 ) + lQdf f) a{q 2 ) 2 ) (13) 



where 



Hereafter we will follow [3] 



4 n/) = 1.9857- 0.115371/ 



4 n/) = 18.244 - 4.216n/ + 0.086n^ 



,/ 2n 2 M<1 2 ) 

d{q ) = -q 



p(q 2 ) = -J 



dq 2 

d(q 2 ) 
Q 2 



dq 2 



(14) 

(15) 
(16) 



Using formula (3) 
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dq 2 da 

~f = ~ (3 a 2 {l + b ia + b 2 a 2 ) 



we get for the function p(q 2 ) the expression 

i I d(a)da 



1 

A> 



1 



d(a)da 



a 2 (l + b\a + b 2 a 2 ) A)&2 ■/ a 2 (a — xi)(a — X2) 
After taking the integral (18) we get 



(17) 



(18) 



4 4 r 1 1 1 

p(a) = —Ina H . < ( h 4c?i + 16d 2 x 1 )ln(a - x ± ) - ( h 4di + 16d 2 x 2 )ln(a - x 2 ) } (19) 



Xi, x 2 are determined by formula (8a). 

The polarization operator is an analytical function with the cut [0,oo]. The polarization operator 
calculated in 3-loop approximation has the physical cut < q 2 < 00 and nonphysical one —A3 < q 2 < 0. 
The contribution of the physical cut in the value R^y+A * s equal to 



R 



QCD 
r,V+A 



= Q\V ud \ 2 S EW [ -4(l-4) 2 (l + 24)(l+r( S )^ + A^ 
J mi mi mi 



phys.cut 



(20) 



where \V u d\ = 0.9735 ±0.0008 is the Cabibbo-Kobayashi-Maskawa matrix element [12], Sew = 1.0194± 
0.040 is the contribution of electroweak corrections [13]. 



f 2 2 

ARW = -24tt 2 ^^ = -0.008. 



mi 



(21) 



is a small correction from the pion pole [3] . The nonphysical cut contribution is equal to 



R 



(QCD) 
r,V+A 



Q\V ud \ 2 S EW ( il(l-4) 2 (l + 24)(l + r( S ))rf £ 
J mi mi mi 



nonphysicalcut 



(22) 



-A 2 



3 Finding of the value A 3 

The value measured in the experiment is 

B{t — > v T + hadrons, S — 0) 



R. 



t,V+A 



GW,J 2 S 



ds 
ml 



„<l\ -\l 'H I - Til 5") 2 

mi 



B(t e u e u T ) 



;i + 2-^)(«i(s) + a^s) + a {s)) - 2^-a (s) 



mi 



mi 



For the value R t ,v+a , the ALEPH [4], OPAL [5] and CLEO [6] Collaborations had obtained 

R T y +A = 3.475 ± 0.022 
Here new value of R Tj s [14,15] is taken into account. 



(23) 



(24) 
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R Tt s = 0.161 ± 0.007 (25) 

The convenient way to calculate the R T in QCD is to transform the integral in the complex s plane 
[16-19] around the circle \s\ = and thus getting a satisfactory agreement with the experiment. 



ds 



R T , V+A = 6m\V ud \*S EW 1 ^(1 - -^) 2 (1 + 2-^)(l + 2-^)n( S ) + Aitf (26) 

J \s\=m^ lil T Iil T Iil T lil T 

There are two values of A 3 , such that R T y + A = 3.475 ± 0.022 (Eq.24). These values are: conventional 
value 



A conv = 61g ± 2 QMeV 



and the alternative new value 

We can calculate R?y+ A also with the hel P of formulae (20,22). At A 3 = 618 ± 29MeV 



A™ w = 1666 ± 7MeV 



R 



QCD 
r,V+A 



= 3.305 ±0.008 



(27) 
(28) 

(29) 



phys.cut 



and 



t>QCD 



= 0.162 ±0.015 (30) 

nonphys.cut 

The sum of integrals on the physical and nonphysical cuts is equal to the integral over the circle (it 
follows from Cauchy theorem) and coincides with the measured value R T y + A (24). The contribution of 
only one physical cut is insufficient to explain the experiment. 
If A 3 = 1666 ± 7MeV 



r>QCD 
^ry+A 



= 3.480 ± 0.0007 



phys.cut 



and 



oQCD 
n r,V+A 



0.0129 ±0.0024 



(31) 



(32) 



phys.cut 



If A 3 = A^ onv -, the nonphysical cut must be taken into account to avoid a discrepancy with the experi- 
ment. If A 3 = A r ^ ew , there is two possibilities. It is possible to omit the contribution of the nonphysical 
cut and to satisfy the requirements of microcausality and unitarity. Alternatively, the contribution of 
the nonphysical cut is taken into account and the requirement of microcausality and unitarity will be 
satisfied only in a future comprehensive theory. 

In refs.[4] R T y and R t a are measured separately. 



R T y = 1.775 ±0.017 



(33) 



R rA = 1.717 ±0.018 

The values R T y and R t> a have been corrected taking into account papers [14,15]. 
In QCD one should have for massless u and d quarks 



(34) 
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Figure 1: Contour integral 



R (QCD) = R (QCD) (35) 

The results of the experiments (33, 34) contradict the formula (35). This contradiction was resolved in 
the paper [2]. 



4 NEW SUM RULES FOR POLARIZATION OPERA- 
TORS. 

To derive the sum rule, let us consider the integral over closed contour from the function W(s)H(s), 
where U(s) is one of the functions U. ( y\s), U^\s), U. { y\s) + IT^ (s), U. ( y\s) - n^s), and W(s) is 
the weight analytical function, which will be chosen later. As a contour, we choose that which contains 
the upper and lower edges of the cut from s\ to S2 and of two circles with radii s\ and S2 (see Fig.l). 
Let us choose the values si = 0.6,0.8, ...2GeV 2 and the values S2 = si + 0.2, si + 0.4, ...3GeV 2 . The 
integral considered through Cauchy theorem is zero. It is does not contain the contribution of power 
corrections 4 and nonphysical cut. As a weight function we choose 

W{s) = {s- Sl )(s 2 -s) (36) 

The sum of integrals over cut edges is 1% J W(s)ImH(s)ds. This sum is equal to the sum of integrals 

with inverse sign over the circles, for which owing to that the weight function vanishes at the points Si 
and S2, one may take Il^ CD \s) instead of the true value n(s). 

Making use of the analytical properties of Il^ CD \s), let us transform the sum of integrals over 
circles into the integral from ImU.^ CD \s) over the cut from s\ to s 2 - Finally, we obtain the following 
sum rule 

S'2 

J W(s)ImU(s)ds 

4 We ignored the a s corrections to the condensates 
contour of integration 



= J W(s)ImU {QCD \s)ds (37) 
. The condensates without a s corrections are the poles off the 
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To compare QCD predictions with experiment, let us introduce the notations 



S2 S2 

U B = J W(s)ImU B (s)ds/ J W(s)ImU ( g CD \s)ds (38) 



Sl Si 

where B = V, A, V + A. The results of the calculations of Ub are given in Table 1. 
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Table 1. 

Comparison of the sum rules (37) with the ALEPH experiment. 
Ub is given by Eq.(38). /mll^s) is obtained from the ALEPH experimental data. ImU B ^ CD \s) is 
calculated by three-loop approximation of QCD at A 3 = 618MeV. si, s 2 are given in GeV 2 . 



S2 


1.2 


1.4 


1.6 


1.8 


2 


2.2 


2.4 


2.6 


2.8 


3 


si = 0.8 


U v 


0.552 


0.489 


0.479 


0.498 


0.534 


0.592 


0.668 


0.745 


0.816 


0.881 


U A 


1.169 


1.398 


1.502 


1.514 


1.465 


1.389 


1.309 


1.229 


1.161 


1.107 


Uv+A 


0.861 


0.942 


0.988 


1.003 


0.995 


0.985 


0.982 


0.981 


0.982 


0.984 


S2 


1.2 


1.4 


1.6 


1.8 


2 


2.2 


2.4 


2.6 


2.8 


3 


Sl = 1 


Uy 


0.438 


0.425 


0.451 


0.494 


0.546 


0.62 


0.71 


0.796 


0.871 


0.939 


U A 


1.491 


1.629 


1.649 


1.594 


1.493 


1.383 


1.279 


1.889 


1.114 


1.059 


Uv+A 


0.967 


1.024 


1.047 


1.039 


1.015 


0.995 


0.998 


0.986 


0.986 


0.988 


«2 


1.4 


1.6 


1.8 


2 


2.2 


2.4 


2.6 


2.8 


3 


si = 1.2 


Uy 


0.429 


0.477 


0.531 


0.592 


0.677 


0.778 


0.869 


0. 945 


1.011 


U A 


1.712 


1.668 


1.563 


1.427 


1.299 


1.188 


1.098 


1.027 


0.978 


U V +A 


1.065 


1.069 


1.042 


1.003 


0.98 


0.976 


0.976 


0.979 


0.983 


S2 


1.6 


1.8 


2 


2.2 


2.4 


2.6 


2.8 


3 


si = 1.4 


Uy 


0.533 


0.586 


0.65 


0.749 


0.862 


0.956 


1.028 


1.091 


UA 


1.601 


1.456 


1.302 


1.172 


1.067 


0.986 


0.927 


0.891 


Uy+A 


1.063 


1.015 


0.968 


0.951 


0.956 


0.963 


0.970 


0.978 


S2 


1.8 


2 


2.2 


2.4 


2.6 


2.8 


3 


si = 1.6 


Uy 


0.642 


0.715 


0.837 


0.962 


1.053 


1.118 


1.173 


UA 


1.299 


1.146 


1.032 


0.944 


0.88 


0.833 


0.816 


Uy+A 


0.964 


0.922 


0.924 


0.944 


0.959 


0.97 


0.98 


S2 


2 


2.2 


2.4. 


2.6 


2.8 


3 


si = 1.8 


Uy 


0.788 


0.956 


1.086 


1.161 


1.209 


1.252 


UA 


1.006 


0.918 


0.847 


0.799 


0.77 


0.765 


U V +A 


0.884 


0.923 


0.958 


0.974 


0.983 


0.922 


S2 


2.2 


2.4 


2.6 


2.8 


3 


si = 2 


Uy 


1.133 


1.217 


1.255 


1.281 


1.312 


U A 


0.833 


0.778 


0.744 


0.728 


0.738 


Uy+A 


0.972 


0.993 


0.996 


0.998 


1.005 



It is seen from Table 1 that QCD predictions agree with experiment for V + A and disagree with 
experiments for V and A separately. The results do not change if one takes the weight function of the 
form (s — Si) n (s 2 — s) n , n = 2, 3, ...10. 



Let us consider V — A. In this case Im{U. 



(QCD), 
V > 



-Uf CD \s)) 



0, while Im(U v (s)-U A (s)) ^ 0. 



Try to eliminate this disagreement with the help of instantons. 



The instanton contribution into Ily^(s) — n^ ; (s) in the model considered in [3] is given by the 
formula (39) [3]: 



r(!) 



n 



a) 

V,inst 



(s) - U% nst (s) = jdpn(p) 



s 



(39) 
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K\ is the Macdonald function. Introduce the notations 

L Exp = J W(s)Im(n$\s)-IL%\s))ds 

Sl 

S2 

L mst = J W(s)Im(n% nst ( s ) - U^ mst (s))ds 

The results of L Exp and L inst calculations are given in Table 2 for n{p) 
1.7GeV-\ n = 1.5 • 1(T 3 GW 4 . 



Table 2. Sl = 0.8GeV 2 . 



s 2 /GeV" 2 


1 


1.2 


1.4 


1.6 


1.8 


LExp 


-0.00042 


-0.022 


-0.107 


-0.286 


-0.554 


Linst 


0.000144 


0.00089 


0.00234 


0.00427 


0.00634 


L Exp/ Li ns t 


-2.907 


-24.147 


-45.8 


-66.82 


-87.34 



s 2 /GeV 2 


2 


2.2 


2.4 


2.6 


2.8 


3 


L>Exp 


-0.878 


-1.198 


-1.442 


-1.572 


-1.563 


-1.402 


L%nst 


0.00815 


0.00933 


0.00955 


0.00853 


0.00619 


0.00209 


L Exp/ Linst 


-107.7 


-128.7 


-151.1 


-184.2 


-256.7 


-670.7 



It is seen from Table 2 that instantons (in the model under consideration) cannot eliminate the 
disagreement between QCD theory and experiment. 

5 NEW QCD PARAMETERS AND ELIMINATION 
OF CONTRADICTIONS 



(40) 
(41) 

= n S(p - p ) Po = 



In my opinion, the only possible way to resolve the discrepancy which follows from the sum rules (37) 
is to change the conventional value A.f nv ~ 600MeV by A^ ew ~ ^OOMeV. Because Si must be larger 
than A3 ~ 2.5GeV 2 , the sum rules (38) become meaningless. At A^ ew ~ lQOOMeV we have a possibility 
to fulfill the requirement of microcausality and unitarity to omit the nonphysical cut. 
The contribution of the physical cut 

R Q T C V D J =0 = S\V ud \ 2 S EW [ ^(l - 4) 2 (1 + + r(s))d8 = 3.483 ± 0.022 (42) 

J mi mi mi 
T T T 

agrees with experiment (24) at 

A 3 = A« ew = (1565 ± 193)MeV (43) 

The value A^ ew in eq.(43) differs from A^ ew in eq.(28) since in (28) we take into account the contribution 
of the nonphysical cut. It is my belief, that the nonphysical cut must be absent. In spite of that the 
nonphysical cut is a consequence of 3-loop in p QCD, if we are able to eliminate this drawback, we must 
do it. But at A 3 ~ QOOMeV we cannot do it while at A 3 ~ lQOOMeV we can. In what follows we put 
A 3 = (1565 ± 193) MeV and omit the nonphysical cut contribuion. 
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6 TRANSITIONS TO A LARGER NUMBER OF 
FLAVOURS 

Let us introduce the notations 



f nf (a) = - + bt ,) Ina- 



^ n/)2 -46. 



(n f ) 
2 



ln{a-x^ !) ) ^ln(a - x { 2 nf) ) 

X2 



X^ 



, n f < 5 (44) 



/n f (a) = - + bi f) 



sjb( n f) 2 -46. 



(«/) 



Zn(a - xJ n/) ) - — j^- /n(x^ n/) - a) 



r W) 2 



, n / = 6 5 (45) 



where 



4 3) = -0.0497 + 0.107i, = )# , 4 4) = -0.0632 + 0.1296i, x { 2 4) = x (A) * , 
xf ] = -0.107 + 0.176i, xf = xf h , xf ] = -0.213, xf = 1.013 
The value a(q 2 ) is found by numerical solution of the equation 

q 2 + iQ\ 



(46) 



A 2 



fn f (a), Q > 0, n f = 3,4,5,6 



(47) 



The function r(s) can be obtained with the help of Eqs. (19,12). 

In the evolution upwards to larger energies the matching of r(q 2 ) at the masses J/ip, T and 2m t is 
performed. 

There are three alternatives: 

1) The nonphysical cut is absent A^ ew = (1565±193)Mel / . The Adler function d(q 2 ) may be written 
in the form 

d(q 2 ) = S 3 \q 2 ) + S 4 \q 2 ) + d^(q 2 ) + d^\q 2 ) (48) 

where 



d (s V) = -9 2 / 



r,{q' 2 )dq' 2 
(q' 2 - q 2 y 



(49) 



is the contribution of the part of the cut with 3 flavours into Adler function. Similarly, 

ml. 

r A (q' 2 )dq' 2 



d {A \q 2 ) = -q 2 J 



(q' 2 - q 2 y 



(50) 



is the contribution of the part of the cut with 4 flavours into Adler function 



3 The sign of the argument in the third logarithm is changed for that to remain on the physical sheet. 
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Am 



2 



« (5 V) = -f I ^% (51) 



(q 12 - q 2 y 

is the contribution of the part of the cut with 5 flavours into Adler function 

r e (q a )dq a 



«<V> - V / gfc^ (52) 



4m? 



is the contribution of the part of the cut with 6 flavours into Adler function. 

The number of flavours for r(q 2 ) on the cut is a certain number in contrast to the number of flavours 
at the point of the complex plane q 2 off the cut. Let us consider q 2 = — m| and find a s (—m 2 z ). 

Return to formula (13). The coefficients d\ and d 2 in (13) are defined for a certain number of flavours. 

Introduce 

d<F\-m 2 z ) = {dfd^\-m 2 z ) + d^S 4 \-m 2 z ) + dfd^\-m 2 z ) + df d^ {-m 2 z )) / d{-m 2 z ) (53) 



dt v \-m z ) = (4 3) rf (3) (-^l) + dfS'\-m\) + dfdP\-m\) + df dP\-m\)) j d{-m\) (54) 
The values d^ n ^{—m z ) have been calculated. 

rf (3) (-m|) = 0.000169, d (4) (-m|) = 0.000823, rf (5) (-m|) = 0.0432 

rf (6) (-m|) = 0.00218, d(-m 2 z )= 0.0464 (55) 

Formula (13) is replaced by 

d(-m%) = 4a(-m|)(l + 4d[ av) (-m 2 z )a(-m 2 z ) + lQd ( 2 av) (-m z )a 2 (-m z )) (56) 
The equation (49) can be solved for a(— m 2 z ) 

a s (-m 2 z ) = 4tt a(-m|) = 0.142 ± 0.004 (57) 
The value a s (m 2 z + iO) is evaluated from (6,8). 

a s {m 2 z + iO) = 0.131 ± 0.003 + (0.037 ± 0.002)i 

\a s (m 2 z + i0)\ = 0.136 ±0.004 (58) 

In a similar way a s (q 2 ) can be calculated at arbitrary q 2 . The values of a s at the interesting points 
are given in Table 3. 
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Table 3. 



The calculation of a s {q 2 ) at different q 2 in approximation of 1-4 loops. The matching of r(q 2 ) at the 
masses of J/ip, T and 2m t is performed. The contribution of nonphysical cut is omitted. 



Approximation 


A 3 /MeV 


AjMeV 


A 5 /MeV 


A 6 /MeV 


« s (0) 


a s (-m 2 T ) 


One loop 


618 ±59 


508 ± 51 


377 ± 40 


191 ±22 


1.396 


0.469 ±0.018 


Two loops 


1192 ± 136 


956 ± 113 


678 ± 86 


301 ±42 


0.895 ±0.001 


0.387 ±0.015 


Three loops 


1565 ± 193 


1257 ± 158 


872 ± 119 


312 ±47 


0.749 ±0.001 


0.379 ±0.013 


Four loops 


1862 ±230 


1503 ± 189 


1064 ± 145 


202 ±31 


0.749 ±0.001 


0.379 ±0.013 



Approximation 


a s (m 2 T + iO) 


a,(-m%) 


a s (ml - + i0) 


One loop 


0.356 ± 0.019 + (0.306 ± 0.017)i 


0.377 ±0.014 


0.276 ± 0.010 + (0.216 ± 0.013)i 


Two loops 


0.355 ± 0.014 ± (0.270 ± 0.014)i 


0.332 ±0.013 


0.265 ± 0.003 + (0.200 ± 0.014)i 


Three loops 


0.365 ± 0.023 + (0.264 ± 0.013)i 


0.322 ±0.009 


0.272 ± 0.015 ± (0.201 ± 0.014)i 


Four loops 


0.363 ± 0.013 + (0.263 ± 0.013)i 


0.326 ±0.012 


0.273 ± 0.016 ± (0.202 ± 0.014)i 



Approximation 


a(ml + + i0) 


a s (-m\) 


a s (my — ± iO) 


One loop 


0.291 ± 0.012 + (0.206 ± 0.013)i 


0.254 ±0.008 


0.206 ± 0.005 + (0.107 ± 006)i 


Two loops 


0.284 ± 0.002 + (0.193 ± 0.014)i 


0.234 ±0.009 


0.192 ±0.004+ (0.100 ±007)i 


Three loops 


0.293 ± 0.017 + (0.199 ± 0.015)i 


0.237 ±0.009 


0.194 ± 0.006 + (0.106 ± 0.008)i 


Four loops 


0.294 ± 0.018 + (0.199 ± 0.015)i 


0.237 ±0.009 


0.194 ± 0.007 + (0.104 ± 0.007)i 



Approximation 


a s (m\ + + i0) 


a s (-m 2 z ) 


a s (ipm 2 + iO) 


One loop 


0.211 ± 0.006 + (0.100 ± 0.005)i 


0.142 ±0.003 


0.141 ± 0.006 + (0.039 ± 0.002)i 


Two loops 


0.199 ± 0.005 + (0.093 ± 0.007)i 


0.137 ±0.002 


0.129 ± 0.003 + (0.036 ± 0.002)i 


Three loops 


0.203 ± 0.007 + (0.098 ± 0.007)i 


0.142 ±0.004 


0.131 ± 0.003 + (0.037 ± 0.002)i 


Four loops 


0.204 ± 0.008 + (0.099 ± 0.007)i 


0.141 ±0.004 


0.131 ± 0.003 + (0.038 ± 0.002)i 



Approximation 


r(mj) 


Ri 


One loop 


0.0462 ±0.0009 


20.856 ±0.017 


Two loops 


0.0458 ±0.0011 


20.848 ±0.021 


Three loops 


0.0464 ±0.0012 


20.860 ± 0.023 


Four loops 


0.0465 ±0.0012 


20.861 ±0.023 



2. The nonphysical cut is taken into account , A^ ew = (1666 ± 7)MeV. In this case the lower 
limit of the integral in eq.(49) is equal to —A3 and in the three- loop approximaion A 4 = 1591Me\^, 
A 5 = 791 MeV, A 6 = 280MeV. Errors in this case are very small. The results of the calculations in the 
three-loop approximation are the following: 

a s (-m 2 T ) = 1.575, a s (m 2 T + iO) = 0.142 + 0.272i, a s {-m\) = 0.484 
a s {ml-0 + i0) =0.179 ± 0.209i a s [m\ + + Oi) = 0.223 + 229i, a s (-m 2 r ) = 0.253 
a s (m 2 r -0 + i0) = 0.202 + O.llOi, a s (m 2 r + + 0i) = 0.189 + 0.093i, a s (-m 2 z ) = 0.139 (59) 
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a s (-m 2 z + i0) = 0.129 + 0.036i, r{m 2 z ) = 0.0457, = 20.845 

3. The nonphysical cut is taken into account, Af nv ' = (618±29)MeU. The results of the calculations 
in 1-4 loop approximaion are given in Table 4. 

Table 4. 



The calculation of a s (q 2 ) at different q 2 in 1-4 loop approximation. The matching of r(q 2 ) at the 
masses of J/ip, T mesons and of 2m T is performed. The contribution of nonphysical cut is taken into 
acccount. 



Approx. 


A 3 /MeV 


A 4 /MeV 


A 5 /MeV 


A 6 /MeV 


a s (-m 2 T ) 


a s (m 2 T + iO) 


One loop 


370 ± 19 


296 ± 16 


211 ± 13 


102 ±7 


0.445 ±0.015 


0.222 + (0.222 ± 007)i 


Two loops 
Three loops 
Four loops 


539 ± 25 
618 ±29 
720 ± 33 


416 ±21 
475 ± 24 
557 ±28 


275 ± 15 
301 ± 16 
359 ± 20 


111 ± 7 
96 ±5 
61 ±4 


0.371 ±0.012 
0.354 ±0.010 
0.375 ±0.012 


0.19 + (0.174 ±0.006)i 

0.186 ± 0.001 + (0.162 ± 0.005)i 

0.184 ± 0.001 + (0.161 ± 0.005)i 



Approx. 


a s (-m%) 


a s (m| - + i0) 


asiml + O + iO) 


One loop 


0.33 ±0.01 


0.212 ± 0.002 + (0157 ± 0.005)i 


0.222 ± 0.002 + (0.148 ± 0.005)i 


Two loops 


0.277 ±0.006 


0.182 ± 0.002 + (0.125 ± 0.004)i 


0.192 ± 0.002 + (0.118 ± 0.004)i 


Three loops 


0.266 ±0.006 


0.177 ± 0.001 + (0.116 ± 0.004)i 


0.187 ± 0.002 + (0.112 ± 0.003)i 


Four loops 


0.273 ±0.006 


0.177 ± 0.001 + (0.116 ± 0.003)i 


0.187 ± 0.002 ± (0.111 ± 0.003)i 



Approx. 


a s (-m\) 


a s (m\ — + iO) 


a s (m\ + + i0) 


One loop 


0.180 ±0.003 


0.180 ± 0.002 + (0.082 ± 0.002)i 


0.184 ± 0.002 + (0.076 ± 0.002)i 


Two loops 


0.189 ±0.003 


0.157 ± 0.002 + (0.066 ± 0.002)i 


0.161 ± 0.002 + (0.061 ± 0.002)i 


Three loops 


0.184 ±0.003 


0.154 ± 0.002 + (0.063 ± 0.002)i 


0.159 ± 0.002 + (0.059 ± 0.001)i 


Four loops 


0.185 ±0.003 


0.154 ± 0.002 ± (0.063 ± 0.002)i 


0.159 ± 0.002 + (0.059 ± 0.002)i 



Approx. 


a s (-m 2 z ) 


a s {m 2 + iO) 


r{m 2 ) 


Ri 


One loop 


0.135 ±0.001 


0.186 ± 0.002 + (0.032 ± 0.001)i 


0.0420 ±0.0004 


20.772 ± 0.008 


Two loops 


0.120 ±0.001 


0.113 ± 0.001 + (0.0274 ± 0.0005)i 


0.0395 ±0.0003 


20.721 ±0.007 


Three loops 


0.118 ±0.001 


0.112 ± 0.001 + (0.0267 ± 0.0004)i 


0.0389 ±0.0003 


20.710 ±0.007 


Four loops 


0.118 ±0.001 


0.111 ± 0.001 + (0.0266 ± 0.004)i 


0.0388 ±0.0003 


20.708 ± 0.007 



Unlike conventional matching procedure at negative q 2 , the matching of r(q 2 ) on the masses of J/^, 
T mesons and of 2m t is performed. The value a s (— m 2 ) is practically independent of the matching 
procedure. I believe, that the first alternative is the best. 



7 Comparison of the calculated values Rt(s) with the mea- 
sured values Re(s). 

The value Rt(s) is calculated by the formula 

R T (s) = 3 £ e 2 q (l + r(s)) (60) 
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The results of the calculations of R{s) in three loop-approximation and of their comparison with ex- 
periments are given in Table 5 for 2 < y/s < A.8GeV and for 12 < y/s < 46.6 GeV in Table 6. The 
calculated values of the function R(s) are in an excellent agreement with the experiment except for the 
resonance region 3.7 < \/s < AA GeV . But the accuracy of measurements of R(s) is insufficient to 
define the value r(s) with a good accuracy. 

Table 5 



Comparison of the calculated values Rt(s) with the measured values Re [21] 



E cm (GeV) 


j } 

JXt 


Ke 


E cm (GeV) 


7") 


Ke 


o r\(~\r\ 

2. (J (J (J 


o on 

2.29 


O 1 O i a r\ 7 i A 1 O 

2.18. ± 0, 07 ± 0.18 


A AO O 

4.033 


O 7A 

3.70 


A OO 1 A 1 77 I A OO 

4.32 ± 0.17 ± 0.22 


2.200 


o oo 

2.28 


O OO 1 A A 7 1 a 1 7 

2.38 ± 0.07 ± 0.17 


A A A A 

4.040 


O 7A 

3.70 


A A A 1 ATT' 1 A 1 A 

4.40 ± 0.17 ± 0.19 


o A a a 

2.400 


O 07 

2.27 


o oo i a a 7 i r\ -\ a 

2.38 ± 0.07 ± 0.14 


A A 7 A 

4.050 


O 7A 

3.70 


A OO 1 n 1 7 I A 00 

4.23 ± 0.17 ± 0.22 


o r a a 

2.500 


O O 7 

2.27 


O OA 1 A AO 1 A 1 7 

2.39 ± 0.08 ± 0.15 


A A/"" 1 A 

4.060 


O 7A 

3.70 


A /** 7 1 A 1 A I A1A 

4.65 ± 0.19 ± 0.19 


O /"'A A 

2.600 


2.26 


O OO 1 A A/" 1 A 1 C 

2.38 ± 0.06 ± 0.15 


4.070 


O 7fl 

3.70 


AAA 1 A OA 1 A 1 A 

4.14 ± 0.20 ± 0.19 


o ta A 

2.700 


2.26 


O OA 1 A A7 1 A 1 O 

2.30 ± 0.07 ± 0.13 


A AO A 

4.080 


O 7fl 

3.70 


/IO/I 1 AOI I A10 

4.24 ± 0.21 ± 0.18 


2.800 


2.25 


O 0^7 1 A A 1 A 1 1 

2.27 ± 0.06 ± 0.14 


A AAA 

4.090 


a 

3.69 


A A 1 A 1 77 I A A O 

4.06 ± 0.17 ± 0.18 


o a a a 

2.900 


2.25 


C\ AO 1 A A 7 1 A 1 O 

2.22 ± 0.07 ± 0.13 


A A AA 

4.100 


O /"* A 

3.69 


O AT 1 A 1 f* 1 A 1 O 

3.97 ± 0.16 ± 0.18 


o a a a 

3.000 


2.25 


O 01 I A A7 I A 1 1 

2.21 ± 0.05 ± 0.11 


A 1 "1 A 

4.110 


O A 

3.69 


O AO 1 A A f* 1 A 1 A 

3.92 ± 0.16 ± 0.19 


o ta a 

3.700 


3.71 


,'\ riO 1 A AO 1 A AO 

2.23 ± 0.08 ± 0.08 


A "1 OA 

4.120 


n a 

3.69 


AAA 1 A O A 1 A OO 

4.11 ± 0.24 ± 0.23 


3.730 


O 71 

3.71 


2.10 ± 0.08 ± 0.14 


A "1 O A 

4.130 


3.69 


O AA 1 A A 7 1 A "1 7 

3.99 ± 0.15 ± 0.17 


3.750 


O 71 

3.71 


O A T 1 A AA 1 A 1 O 

2.47 ± 0.09 ± 0.12 


A 1 A A 

4.140 


O A 

3.69 


O OO 1 A A 7 1 A 1 O 

3.83 ± 0.15 ± 0.18 


3.760 


3.71 


O 77 1 A 1 1 I A 1 O 

2.77 ±0.11 ± 0.13 


,4 1 7 A 

4.150 


a 

3.69 


A O A 1 A 1 O 1 A1A 

4.21 ± 0.18 ± 0.19 


O ^ O A 

3.764 


3.71 


O OA 1 A OT 1 A OA 

3.29 ± 0.27 ± 0.29 


A 1 a 

4.160 


O /"* A 

3.69 


A A O 1 A 1 7 1 A 1 O 

4.12 ± 0.15 ± 0.16 


o ^ o o 

3.768 


O ^71 

3.71 


O OA 1 A OO 1 A 07 

3.80 ± 0.33 ± 0.25 


A A 7A 

4.170 


O /"* A 

3.69 


A 1 O 1 A 1 7 1 A1A 

4.12 ± 0.15 ± 0.19 


O 77 A 

3.770 


O 7 1 

3.71 


077 1 A 1 A I A1A 

3.55 ± 0.14 ± 0.19 


/I 1 OA 

4.180 


a 

3.69 


A AO 1 A 1 7 1 AIO 

4.18 ± 0.17 ± 0.18 


3.772 


3. 71 


OIO I A O A I A OO 

3.12 ± 0.24 ± 0.23 


/I 1 AA 

4.190 


a 

3.69 


A r\A 1 r\ a a 1 r\ a a 

4.01 ± 0.14 ± 0.14 


3.776 


O 71 

3. 71 


O O t • I A O ( • 1 A 1 A 

3.26 ± 0.26 ± 0.19 


A OA A 

4.200 


a 

3.69 


OOT 1 A 1 Z 1 I A 1 r* 

3.87 ± 0.16 ± 0.16 


3.780 


O ^71 

3.71 


O OO 1 A 1 O I A 1 O 

3.28 ± 0.12 ± 0.12 


A OI A 

4.210 


O /"* A 

3.69 


O OA 1 A A O 1 AIT 

3.20 ± 0.16 ± 0.17 


o 7 a a 

3.790 


3.71 


O ( ' O 1 ATI 1 A 1 A 

2.62 ± 0.11 ± 0.10 


/I OO A 

4.220 


O A 

3.69 


O ( ' O 1 A A 7 1 A OA 

3.62 ± 0.15 ± 0.20 


3.810 


3. ( 1 


Z.08 ± 0.10 ± U.lz 


a oon 
4.ZO0 


3.09 


01 _i_ n 1 _i_ n 1 r 
6.11 ± 0.13 ± 0.15 


3.850 


3.70 


2.47 ±0.11 ±0.13 


4.240 


3.69 


3.24 ±0.12 ±0.15 


3.890 


3.70 


2.64 ±0.11 ±0.15 


4.245 


3.69 


2.97 ±0.11 ±0.14 


3.930 


3.70 


3.18 ±0.14 ±0.17 


4.250 


3.69 


2.71 ±0.12 ±0.13 


3.940 


3.70 


2.94 ±0.13 ±0.19 


4.255 


3.69 


2.88 ±0.11 ±0.14 


3.950 


3.70 


2.97 ±0.13 ±0.17 


4.260 


3.69 


2.97 ±0.11 ±0.14 


3.960 


3.70 


2.79 ±0.12 ±0.17 


4.265 


3.69 


3.04 ±0.13 ±0.14 


3.970 


3.70 


3.29 ±0.13 ±0.13 


4.270 


3.69 


3.26 ±0.12 ±0.17 


3.980 


3.70 


3.13 ±0.14 ±0.16 


4.280 


3.69 


3.08 ±0.12 ±0.15 


3.990 


3.70 


3.06 ±0.15 ±0.18 


4.300 


3.69 


3.11 ±0.12 ±0.12 


4.000 


3.70 


3.16 ±0.14 ±0.15 


4.320 


3.69 


2.96 ±0.12 ±0.14 


4.010 


3.70 


3.53 ±0.16 ±0.20 


4.340 


3.69 


3.27± 0.15 ±0.18 


4.020 


3.70 


4.43 ±0.16 ±0.21 


4.350 


3.69 


3.49 ±0.14 ±0.14 


4.027 


3.70 


4.58 ±0.18 ±0.21 


4.360 


3.68 


3.47 ±0.13 ±0.18 


4.030 


3.70 


4.58 ±0.20 ±0.23 


4.380 


3.68 


3.50 ±0.15 ±0.17 
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Table 5. Continuation 



E cm (GeV) 


Rt 


Re 


4.390 


3.68 


3.48 ±0.16 ±0.16 


4.400 


3.68 


3.91 ±0.16 ±0.19 


4.410 


3.68 


3.79 ±0.15 ±0.20 


4.420 


3.68 


3.68 ±0.14 ±0.17 


4.430 


3.68 


4.02 ±0.16 ±0.20 


4.440 


3.68 


3.85 ±0.17± 0.17 


4.450 


3.68 


3.75 ±0.15 ±0.17 


4.460 


3.68 


3.66 ±0.17 ±0.16 


4.480 


3.68 


3.54 ±0.17 ±0.18 


4.500 


3.68 


3.49 ±0.14 ±0.15 


4.520 


3.68 


3.25 ±0.13 ±0.15 


4.540 


3.68 


3.23 ±0.14 ±0.18 


4.560 


3.68 


3.62 ±0.13 ±0.16 


4.60 


3.68 


3.31 ±0.11 ±0.16 


4.80 


3.67 


3.66 ±0.14 ±0.19 



Table 6. 

Comparison of the calculated values Rt{s) with the measured values R E [22-24] 



s/Gev 


Rt 


Re 


y^/GeV 










[22] 








14.0 


3.92 


4.10 ±0.11 ±0.11 


29.93 


3.88 


3.55 ±0.40 ±0.11 


22.0 


3.89 


3.86 ±0.12 ±0.11 


30.38 


3.87 


3.85 ±0.19 ±0.12 


33.8 


3.87 


3.74 ±0.10 ±0.10 


31.29 


3.87 


3.83 ±0.28 ±0.11 


38.3 


3.86 


3.89 ±0.10 ±0.09 


33.89 


3.87 


4.16 ±0.10 ±0.12 


41.5 


3.86 


4.03 ±0.17 ±0.10 


34.50 


3.87 


3.93 ±0.20 ±0.12 


43.5 


3.86 


3.97 ±0.08 ±0.09 


35.01 


3.87 


3.93 ±0.10 ±0.12 


44.2 


3.86 


4.01 ±0.10 ±0.08 


34.45 


3.87 


3.93 ±0.18 ±0.12 


46.0 


3.86 


4.09 ±0.21 ±0.10 


36.38 


3.87 


3.71 ±0.21 ±0.11 


46.6 


3.86 


4.20 ±0.36 ±0.10 


40.32 


3.86 


4.05 ±0.19 ±0.14 






[23] 


41.18 


3.86 


4.21 ±0.22 ±0.14 


29 


3.88 


3.96 ±0.09 


42.55 


3.86 


4.20 ±0.22 ±0.14 






[24] 


43.53 


3.86 


4.00 ±0.20 ±0.14 


12 


3.93 


3.45 ±0.27 ±0.13 


44.41 


3.86 


3.98 ±0.20 ±0.14 


14.04 


3.92 


3.94 ±0.14 ±0.14 


45.59 


3.86 


4.40 ±0.22 ±0.15 


22 


3.89 


4.11 ±0.13 ±0.12 


46.47 


3.86 


4.04 ±0.24 ±0.14 


25.01 


3.88 


4.24 ±0.29 ±0.13 








27.66 


3.88 


3.85 ±0.48 ±0.12 
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8 COMPARISON WITH a s (s) OBTAINED FROM THE 
SUM RULES 

Perturbative corrections to two measurements, namely, the Gross-Llewellyn Smith sum rule [7] for the 
deep inelastic neutrino scattering and the Bjorken sum rule [8] for polarized structure functions, have 
been determined: 

a s (-3GeV 2 ) = 0.28 ± 0.035(stat) ± 0.050(sys.) [12, 25] (61) 

and 

a s (-2.5GeV 2 ) = 0.375lg;{jjj? [12, 26 - 28] (62) 
The results of the calculations give 

a s (-3GeV 2 ) = 0.381 ± 0.013 (63) 
a s (-2.5GeV 2 ) = 0.39 ± 0.013 (64) 

9 Calculation of Ri 

The value Ri = Y[Z — > hadrons) /T(Z — > leptons) is parametrized by the latest version of ZFITTER 
[29] 

Ri = 19.934( 1 + 1.045(— ) + 0.94(— ) 2 - 15(— ) 3 ) (65) 

y 7T 7T 7T J 

This is the conventional method. In this parametrization 

M H = 300Gel/, M t = HA.lGeV, 0.1 < a s < 0.13. 
In the method under investigation Eq.(65 ) should be changed by the formula 

Ri = 19.934(1 + r(m|)) (66) 

The calculated values of r{m 2 ) and of Ri are presented in Tables 3,4 in 1-4 loop approximation. The 
value Ri can be compared with the measurement [12]. 

T(Z -> hadrons) = (1744.4 ± 2) MeV (67) 
T(Z -> /eptons) = (83.984 ± 0.086) Me^ (68) 

and 

Ri = 20.771 ± 0.045 (69) 
The value Ri (69) does not contradict Ri from Tables 3,4. 
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10 On analiticity of a s (q 2 ) 



In the interesting papers [30] the renormgroup was combined with analiticity of a s (q 2 ). It was 
assumed that a s {q 2 ) is an analytic function of q 2 in the whole complex q 2 plane with a cut along the 
poitive q 2 semiaxis. In particular, it was obtained that a s (0) is universal and independent of the number 
of loops, A 3 and 

a s (0) = 4vr//3j 3) = 1.396 (70) 

The analogous formula had been also obtained in paper [1]. As is seen from Table 3, eq.(70) is valid 
only in one-loop approximation. 

If a s (q 2 ) has correct analiticity, the Adler function will have correct analiticity too, but a reverse 
statement is invalid. If Adler function has correct analiticity, function a(q 2 ) may have, generally speak- 
ing, additional singularities. This statement will be evident, if one considers eq.(13) as an equation 
relative to a(q 2 ). If only eq.(13) is solved by expansion in a(q 2 ), analiticity a(q 2 ) will be the same as 
the Adler function. However, using of d(q 2 ) expansion in a(q 2 ) at small q 2 seems to be doubtful. 

11 Conclusion. 

Concluding, let us formulate the main results of this paper. 

1) There are two and only two values of A3, at which R T y + A = 3.475 ± 0.022, one conventional 
value Af nv = (618 ± 29) MeV (A 3 is defined by eq.(7)) and the other, found in this paper, A^ ew = 
(1666 ± 7) MeV. 

2. The renormgroup calculation leads to appearance of a nonphysical cut in the Adler function. 
In complete theory, where everything is taken into account, the nonphysical cut must be absent. The 
question arises, if it is possible to neglect nonphysical cut at the present situation with the theory. At 
conventional value A 3 = Af nv it is impossible. The new value, A 3 = A^ ew is more preferable than 
A 3 = Af nv , since at A 3 = A 3 ew the contribution of nonphysical cut into R T y + A is practically absent. 

3) At A 3 = Af nv there is an essential disagreement between the ALEPH experiment and new 
obtained sum rules, which follow only from analytical properties of the polarization operator. This 
disagreement disappears if A 3 = A^ ew . 

4) In 1-4 loop approximation all calulations are made exactly, without rr/lnjj expansion. 

5) At A 3 e-lu = (1565 ± 193) MeV the nonphysical cut may be omitted, the polarization operator has 
correct analytical properties and R T y + A = 3.475 ± 0.022. But in this case QCD parameters essentially 
differ from conventional values. 
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